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Abstract. Pakistan has faced numerous natural disasters like floods, earthquakes, landslides 
and environmental degradation which severely affects the Pakistan’s economy and results in 
various problems like causalities, diseases, water stress and severe damages (e.g., houses, public 
infrastructure and agricultural land erosion). There is a lack of systematic approaches to 
analyze pre and post damage assessment for estimation of exact loses and the total cost for 
rehabilitation of damaged infrastructure in an efficient way. There exist a variety of 
mechanisms but GIS based flood mapping is considered the most efficient to manage the 
flood situation. This study is focused on evaluation of flood affected areas especially in Punjab 
using WorldView-2, 8-band multi-spectral imagery by applying Remote Sensing (RS) and GIS 
techniques. The research area is comprised of Kot Addu and Muzaffargarh Districts in Punjab 
province of Pakistan that faced a catastrophic super flood of 2010. The WorldView-2, Quick 
Bird and multispectral satellite imagery are capable of making better decisions and assessment 
of flood effected area accurately. RS and GIS techniques can achieve the objectives and 
significant analyses through visual interpretations. These techniques are also used to identify 
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the flood affected regions. The study site was examined by applying supervised classification 
on the basis of the training areas which were obtained during the field surveys in the study 
site. Supervised classification determines that 16900.96 Hectors of agriculture land was 
damaged while Sparse Riverine Forest had the area 44.52 hectors. The damaged built-up area 
was 1805.78 Hectors. RS and GIS techniques are efficient for flood mapping. 


Keywords: Multi-spectral bands, GIS, RS, Visual Interpretation Elements, Supervised 
Classification, Worldview-2 and Quick Bird. 


Introduction. 


Natural disasters such as flood, tsunami, typhoon forest fire, earthquake, volcanic 
eruption and debris flow often cause interference to cultural environment and physical 
structures, houses and infrastructure [1]. Flooding creates a marvelous influence on human 
life and has great worldwide importance regarding hydrological cycle [2]. Each year natural 
hazards, such as earthquake, typhoon, flood, and debris flow etc., cause tremendous loss of 
properties and human lives in the world [3]. The hazardous floods are the most shocking, 
widespread and frequent, amongst other natural disasters [4,5]. Flood vulnerability represents 
the potential of flood that portends floodplain inhabitant’s lives and their assets and 
agricultural land [6]. Floodplain is relatively plane and low-lying area where water flow over 
the river banks [7,8,9,10]. As floods are the natural phenomenon but flood hazards cause to 
devastate the social life of inhabitants, business, housing, and agricultural activities, which 
often severely damaged during flood [6,11,12,15]. 

In Pakistan, floods are classified on the basis of pivotal factors e.g., riverine flood, dam 
break flood and storm surge floods [14]. Riverine flood happens due to the inflow of water to 
a river, when water surpasses to river’s channel capacity [12]. In Pakistan, flood history from 
1950-2015, highlights many high medium and low-level floods along Indus River and its 
connecting channels but the flood event occurred in 2010 was super flood which resulted in 
major damages to adjacent districts along Indus and Chenab River. According to federal flood 
commission of Pakistan about 1,985 people lost their lives, 1.6 million houses were smashed, 
17,553 villages and a total area of 160,000 km^ were affected by the flood occurred in 2010. 
This flood damaged infrastructure, roads, irrigation and drainage system throughout the 
country and effected the social life of about 20 million people across the country. 

In some scenario, flood safety works fascinate high value land use in flood hazards 
areas but provide false sense of security to floodplain dwellers and their properties therein 
because these works often get breached during high flood flow [15, 12]. The main factor of 
breaks in river channels is the sand molding which is due to shocking floods [16]. 

Pakistan is struggling to manage natural hazards for survival from nastiest catastrophic 
conditions which may occur in future. About one fifth of Pakistan’s area was affected by flood 
water in 2010 which pushed thousands of families to displace for their survival and leaving 
millions as homeless. The low-lying Indus belt in the Punjab Pakistan represents a flood prone 
zone. This flood hit the areas along Chenab and Sutlej River. It destroyed the existing 
infrastructure completely. Both the cultural and economic standards of life had been drained 
out in rain. Flood 2010 covered the complete study site in comparison to floods of 2014, 2015 
as shown in figure 1. 
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Figure 1: Flood Extent 2010, 2014 and 2015 in the study site. 
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The management of a disaster in efficient ways and the rehabilitation activities are 

significant for the prosperity of affected people. Intime management of Hazards and damage 
assessments are key indicators for defining actions that communities can take to reduce socio- 
economic losses [3]. Adequate understating of the flood assessment can support the decision 
makers in choice of engineering or non-engineering measures for flood damages [17]. 
There may be different measures that can be taken for hazard management in a region, such 
as land use planning and management, construction of engineering structures and flood 
monitoring ot warning systems, etc. The effectiveness of these alternatives can be more 
properly evaluated if the compatisons are based on the risk reductions from the 
implementation of that measure [18]. 

In order to estimate the damages occurred due to flood hazards, optical data is widely 
used [19,20,21,22,23,24,25,26]. For land cover identification and assessment accurately, 
satelite imagery is considered as dominant tool [27]. The satellite imagery having spatial 
resolution more than 5 meter creates great impact to estimate the accurate extent of flood [28]. 
The spatial resolution provides great value to understand the variations in topography and 
various natural and manmade activities [29]. In order to execute the classification for land 
cover mapping, usually multispectral imagery is used [30]. For image exploration, classification 
is considered as the most significant part [31] used to retrieve the information about different 
land cover classes from the image. To take trainee areas for supervised classification, different 
band combinations and interpretation keys e.g., texture, tone pattern, shadow, 
lightness/darkness, associate and situation of land cover features in image ate widely used [30]. 

A variety of classification techniques are available for Landuse Landcover 
reorganization [32]. Supervised classification is one of these techniques widely used in various 
applications. Supervised classification is based on detection algorithm using pixels from 
known reference samples, usually located within a scene, as a key for comparison to other 
pixels in the same scene [33]. 

Assessment of flood affected areas was examined using WorldView-2 8-band sample 
imagery. The Hyper-spectral imagery has capacity to make better decision and assessment of 
the focused area. GIS and RS techniques are used in this study to achieve the objectives and 
relevant analysis. Current status of the Land Cover and Land Use (LULC) is estimated by 
using supervised classification on the basis of the training areas which were selected during 
the field survey through GPS. 

LULC classification can be seen on a continuum, starting with a basic estimation of 
land cover through broad categories, like farmland and urban areas, to feature extraction, like 
road networks, buildings, and trees [29]. The visual interpretation technique is the most precise 
approach [34] and is used for comparison of supervised classification results and identification 
of flood affected areas. WorldView-2 imagery has specific cartelistic with 1.86 m multispectral 
resolution by using False Color Composite (FFC), panchromatic 0.46m imageries. The data 
was compared with the Quick Bird pre-imagery and digital photographs to sort out the damage 
areas. 

Digital Globe World View-2 has capacity to discriminate features in the study site in a 
better way through bathymetric analysis and vegetative analysis with more accuracy as it is 
possessing 8-spectral bands with additional four bands. With hyper spectral and coastal band 
to asses and estimate the flood extent, damage, harm etc. World View-2 consists of one 
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panchromatic and two multispectral imaging sensors; MS1 includes Blue, Green, Red, and 
Near-IRI and MS2 includes coastal Blue, Yellow, Red Edge and Near-IRI [35]. 


The objective of this study was to utilize Remote Sensing (RS) techniques with the 
available Geographic Information System (GIS) to develop a rehabilitation model consisting 
of transparent information regarding total flood affected area & identification of demolished 
household, damaged infrastructure, flooded crops to achieve the level of better decision 
making to help the responsible authorities for developing, designing, and operating flood 
control infrastructure and preparing aid and relief operations for high-risk areas during future 
floods. 

Material and Methods. 

Investigation site. 

The study was conducted in Tehsil Kot Addu of District Muzaffargarh, Punjab. The 
spatial location of study area is 70.929°E & 30.944°N. District Muzaffargarh is among the 
severely damaged districts during the flood because it is surrounded by two rivers (Chenab 
and Indus). Both rivers were fully flooded which dragged most of the district area under water 
as shown in Figure 2. 

Kot Addu is the severely damaged tehsil of the district Muzaffargarh as shown in the 
table 1. The number of houses damaged in Kot Addu reached up to 71,104 along with six 
fully affected basic health units leading to miseries and displacements [18]. 


Table 1: Summary of Total Affected Area of District Muzaffargarh 


Name of Total Total Affected Several Affected Moderately Affected 
Tehsil UCs Villages Villages UCs (>50%) UCs (<50%) 
Muzaffargarh 35 426 317 12 17 
Kot Addu 28 341 174 21 2 
Alipur 14 124 E pP 7 
Jatoi 16 93 35 0 5 
Total 93 984 598 37 31 
Data used 


For this study hyperspectral imagery of spatial resolution 1.86m and panchromatic 
of spatial resolution 0.46m of World View-2, Quick Bird Scene and the GIS data (Source : 
The Urban Unit) are used. About 400 Square kilometer area (38913.87 Hectors) of Tehsil Kot 
Addu, district Muzaffargarh was acquired from Digital globe for the research purpose. World 
view- 2 has 8 bands with 4 extra bands. The spectral range of the 8-bands is given as in Table- 
2. 

Table 2: WorldView-2 Multispectral Range 


d Band Name Spectral Range 
1 Coastal Blue 400 nm-450 nm 
2 Blue 450 nm-510 nm 
3 Green 510 nm-580 nm 
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4 Yellow 585 nm-625 nm 
5 Red 630 nm-690 nm 
6 Red-Edge 705 nm-745 nm 
7 NIR1 770 nm-895 nm 
8 NIR2 860 nm-1040 nm 
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Figure 2: Location Map of Study Area 
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Methodology 

The data was provided by the Digital Globe in patches. The mosaicking process was 
performed in Erdas Imagine to form a single scene of the WorldView-2 imageries. Moreover, 
the panchromatic image was diffused with multispectral imagery to get more clarity and high 
quality for visual interpretation. To identify the damaged features, visual interpretation 
technique is performed by using FCC on the multispectral 8-band imagery of spatial resolution 
1.86m and compare with WorldView-2 panchromatic imagery of spatial resolution 0.46m, 
Quick Bird scenes of spatial resolution 0.6m along with digital photographs. About 170 
categories of different LULC were identified using supervised classification on the basis of 
training areas and further recoded into 13 broad categories. GPS coordinates and the digital 
photographs of these damaged roads, drainage network, abadies, education institutes, health 
units, bridges, different kinds of agricultural and water areas were obtained during the field 
survey for verification of land use classification. The methodology flow chart is given in figure 
3. The thematic map of training/samples points with WorldView-2 multispectral imagery is 
shown in the form of figure 4. 
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Figure 3: An Overview of the Methodology 
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Figure 4: Training/Sample Points Map on the WorldView-2 Multispectral Image in FCC 
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Results and Discussion 

Supervised classification was executed with the help of the sample points which were 
collected during the field survey. By using coastal blue band, the deep water and the 
sedimented water 1s identified. Yellow and Red-Edge spectral bands have the major role in 
identification of different kinds of the agriculture field. Even by using NIR2 spectral band, the 
damaged agriculture field is identified which has the dominant area than the other classes. 
Riverine forest, sparse forest, plantation/orchards, land soil and built-up areas are also 
identified in the classification results. About 13 classes during supervised classification results 
are shown in the thematic map as in figure 5. The area wise detail of the 13 classes is shown 
in the table 3. 

Table 3: Supervised Classification Result 


Sr. No. Class Name Area (Hector) 
1 Healthy Agriculture Fields 540.69 
2 Harvested Fields 974.89 
3 Damage Agriculture Field 16900.96 
4 Ploughed Fields 3684.62 
5 Plantation/Orchards 455.03 
6 Sparse Riverine Forest 44,52 
7 Dense Riverine Forest 190.34 
8 River Bed 2385.81 
9 Land Soil 2070.64 
10 Wet Land 7330.79 
11 Deep Water 440.79 
12 Sedimented Water 2088.99 
13 Built-up 1805.78 

Total 38913.87 


The result show that 16900.96 Hectors area out of 37108 Hectors was damaged agriculture 
land which is highest than the harvested field, healthy agriculture and ploughed field, while 
Sparse Riverine Forest has the area 44.52 hectors. 


6.2: Visual interpretation 


Visual interpretation technique is used to identify the damaged features using FCC on the 
multispectral 8-band imagery of spatial resolution 1.86m which was compared with 
WorldView-2 panchromatic imagery of spatial resolution 0.46m, Quick Bird scenes of spatial 
resolution 0.6m. 


The direct tangible damages and intangible damages were observed by the flood 2010. 
Direct tangible damages include the physical damage of the houses and infrastructures, while 
intangible damages include the health impacts; physical and mental impacts. The different 
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tangible damaged features; school, village, road, agriculture, bridge, and irrigation channel are 
shown in figures 6,7, 8, 9, 10, 11 and 12. The school is worst damaged during the flood 2010 
as shown in figure 7 (a, b, c and d). Village is highlighted in the figure 8 (a, b, c and d), which 
is also badly affected during the flood 2010. The infrastructure was consist of; road, bridge 
and irrigation channel are also badly affected during the flood 2010. The significance of the 
flood can be seen in such kind of visual interpretations. The crops were badly affected in flood 
2010 which is shown in the figure 10 (a, b, c, and d). The standing flooded water can be seen 
in digital photograph. 
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Figure 5: Supervised Classification Map 
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Figure 6: Supervised Classification Chart 


In case of intangible damages through flood 2010, only 2 persons died and 35 persons were 
injured while 2127 cattle head lost in district Muzaffargarh according to PDMA reports 
[26,37,38,39]. 


Rehabilitation 


In order to develop rehabilitation model for the flood affected infrastructures, damage 
houses and crops, the provincial government and federal government took initiative and 
compensated the losses of the flood affected people. The provincial government repaired the 
flood affected infrastructures (road, school, irrigation channels, bridges, masjids, health 
clinics/dispensaty, tube wells, etc) and prepared the model villages for the flood affected 
people as shown in figure 13. The houses in the model villages were allotted after the 
verification to the flood affected people. The school, dispensaries, tevta centers, livestock 
shed/veterinary dispensary, park and community center are also included. 


Now National Disaster Management Authority (NDMA), Provincial Disaster Management 
Authority (PDMA), Irrigation Department, Planning & Development department and district 
level government are planning to maintain the disaster response plans, risk vulnerability 
assessment and preparing the spatial decision support system (SDSS) to mitigate the 
rehabilitation activities before future flood. 
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Figure 7: Damaged GGPS, Shauhrat wala UC 1 School Building is shown temporally in a, 
b, c and d 


(c) Pictorial View 


(a) World View-2, Multispectral (b) World View-2, Panchromatic 
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Figure 8: Damaged Basti Korai Wala is shown temporally in a, b, c and d and d 
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(a) World View-2, Multispectral (b) World View-2, Panchromatic 
Band Combination: 8, 6, 2 Date Acquisition: 20-12-2010 


Date Acanicitinn: 20-12-2010 


(c) Digital Photograph (d) Quick Bird, 0.6m 
Date Acquisition: 21-10-2009 


Figure 9: Damaged Taunsa Road is shown temporally in a, b, c and d 


(b) World View-2, Multispectral 
Band Combination: 8, 6, 2 
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(c) Digital Photograph ) Quick Bird, 0.6m 
Date Acquisition: 21-10-2009 


Figure 10: Damaged Agricultural at Haala village is shown temporally in a, b, c and d 
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(c) World View-2, Multispectral (b) World View-2, Panchromatic 
Band Combination: 8, 6, 2 Date Acquisition: 20-12-2010 
Date Acquisition: 20-12-2010 


(d) Quick Bird, 0.6m 
(c) Digital Photograph Date Acquisition: 21-10-2009 


Figure 11: Damaged Bridge is shown temporally in a, b, c and d 


(d) World View-2, Multispectral (b) World View-2, Panchromatic 
Band Combination: 8, 6, 2 Date Acquisition: 20-12-2010 


a - 


(d) Quick Bird, 0.6m 


(c) Digital Photograph 
Date Acquisition: 21-10-2009 


Figure 12: Damaged Irrigation channel is shown temporally in a, b, cand d 
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Figure 13: Location Map of Model Villages in Study Area 
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Conclusions 

For disaster rehabilitation plan, it is of great importance to obtain classification results with 
precision. Worldview-2 multispectral image has the great potential to identify the existing 
landuse features by using the extra 4 bands besides its high spatial resolution. 'This imagety 
can be used for vatious application and obtain a fruitful result by using this study for different 
application related to natural disaster estimation, assessment, management, and 
redevelopment of affected areas. The combination of 8, 6, and 2 was found to be useful for 
flood related application. It made the differentiation of vegetation gone under water very clear. 

The red edge band and NIR2 enhanced the vegetation even under water and on the other 

hand blue band helped to identify water in fields. The visual interpretation of the affected area 

also worked well. The settlements and agricultural lands were efficiently mapped. The 
implementation of coastal blue band helped to identify the sedimented water in comparison 
to the deep blue water during supervised classification. 

Limitation of the study 

The study focuses on the damage assessment of the flood affected area in 31“ July 2010 to 

27^ August 2010. To properly analyze the pre-flood and post flood situations of the area, a 

temporal imagery from world view-2 was required. Due to non-availability of the pre-flood 

imagety, an imagery of 0.6 m from the quick bird is used for temporal analysis. 

Further Implications 

8-band imagery has good capacity to identify the damage features, vegetation analysis, and 
water analysis by using extra 4 bands. So, disaster management specifically flood can be 
efficiently mapped on al large scale, e.g., whole Punjab province after flood. 
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